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We have found an exact phase diagram of the Potts model with next nearest neighbor interactions 
on the Bethe lattice of order two. The diagram consists of five phases: ferromagnetic, paramagnetic, 
modulated, antiphase and paramodulated, all meeting at the Lifshitz point i.e. p — 1/3. We report 
on a new phase which we denote as paramodulated, found at low temperatures and characterized by 
2-periodic points of an one dimensional dynamical system lying inside the modulated phase. Such 
a phase, inherent in the Potts model has no analogues in the Ising setting. 
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The (/-state Potts model is one of the most studied 
models in statistical mechanics due to its wide theoreti- 
cal interest and practical applications 
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Potts model @] was introduced as a generalization of the 
Ising model to more than two components and presently 
encompasses a number of problems in statistical physics 
(see, e.g. Q). The model is structured richly enough 
to illustrate almost every conceivable nuance of the sub- 
ject. While some exact results regarding certain proper- 
ties of the model are known, most of them are based on 
approximation methods. In fact, to date, no analytical 
solutions on standard lattices are known to exist. Inter- 
estingly enough however, investigations on phase transi- 
tions of spin models on hierarchical lattices provides ex- 
act calculations of various physical quantities 0, 0, H, Q ■ 
Such studies on the hierarchical lattices begun with de- 
velopment of the Migdal-Kadanoff renormalization group 
method where the lattices emerged as approximants of 
the ordinary crystal ones. It is believed that several 
among its interesting thermal properties could persist for 
regular lattices, for which the exact calculation is by far 
intractable. In [H,[n|, the phase diagrams of the g-state 
Potts models on the Bethe lattices were studied and the 
pure phases of the the ferromagnetic Potts model were 
found. The Bethe lattices were fruitfully used, providing 
a deeper insight into the behavior of the Potts models. 

On the other hand, compared to the Ising models 
with competing interactions the Potts models with 
such interactions on regular and trees are less studied 
P, [H, Hi HJ. In [Hi, the phase diagram for the q- 
state Potts model is constructed by means of the low- 
temperature expansion technique. An infinite set of 
phases appears with the bifurcating structure resembling 
the complete Devil's staircase. In [23|, [24 [ the three-state 
Potts model with antiferromagnetic nearest-neighbor and 
ferromagnetic next-nearest-neighbor interaction was in- 
vestigated within a mean-field theory. 

To the best knowledge of the authors, g-state Potts 
model with competing interactions on the Bethe lattices 
are not well studied [2a] . In this Letter we present a 
phase diagram of the three-state Potts model with next 



nearest neighbor interactions on a Bethe lattice of order 
two. The similarity of results obtained for models de- 
fined on Bethe lattices and on crystal lattices presents 
a strong motivation for the study of models on Bethe 
lattices, since statistical mechanics on such lattices pre- 
sets many simplifying aspects that are absent in models 
defined on crystal lattices. One of the useful ways to in- 
vestigate models defined on trees is to formulate them as 
dissipative mapping problems which allows us to use the 
techniques of the theory of dynamical systems. For the 
models defined on crystal lattices, such a method does 
not lend itself to a simpler solution of the problem since 
the metastable confi gur ations correspond to unstable or- 
bits of the mapping [28| . 

Recall that the Bethe lattice T k of order k > 1 is an 
infinite tree, i.e., a graph without cycles with exactly k+1 
edges issuing from each vertex. Let T k = (V, A), where 
V is the set of vertices of T k , A is the set of edges of T k . 
Two vertices x and y are called nearest neighbors if there 
exists an edge I £ A connecting them, which is denoted 
by / =< x, y >. The distance d(x,y),x,y £ V, on the 
Bethe lattice, is the number of edges in the shortest path 
from x to y. For a fixed x° £ V we set V n — {x £ 
V\d(x,x°) < n} and L n denotes the set of edges in V n . 
For the sake of simplicity we put \x\ — d(x 1 x°), x £ V . 
Two vertices x,y £ V are called the second neighbors if 
d(x,y) = 2. The second neighbor vertices x and y are 
called next nearest neighbor if \x\ ^ \y\ and is denoted by 
> x,y <. 

In this letter, we will consider a semi-infinite Bethe lat- 
tice r?_ of order 2, i.e. an infinite graph without cycles 
with 3 edges issuing from each vertex except for a; which 
has only 2 edges. Considering the three-state Potts model 
with spin values in $ — {1,2,3}, the relevant Hamilto- 
nian with next nearest neighbor interactions has the form 



<x,y> 



(1) 



>x,y< 



where J p , J\ £ R are coupling constants and 5 is the 
Kronecker symbol. In what follows, we consider the case 
where J\ > and J p < 0. 
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In order to produce the recurrent equations, we con- 
sider the relation of the partition function on V n to 
the partition function on subsets of V n -i- Given the 
initial conditions on V\, the recurrence equations indi- 
cate how their influence propagates down the tree. Let 
Z^(ii,io,i2) be the partition function on V n where the 
spin in the root x° is %q and the two spins in the pro- 
ceeding ones are i\ and i%, respectively. There are 27 
different partition functions Z^ 1 ' (i\,i§,%<2) and the par- 
tition function Z^ n ' in volume V n can the be written as 
follows 



Z (n) = £ Z^(h,i Q ,i 2 ). 
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as follows 
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As shown in 27] one can select only five indepen- 
dent variables (1,1,1), Z ( - n \2,l,2), ZW(1,2,1), 
Z(™)(2,2,2), Z(™)(3,2,3) and with the introduction of 
new variables 
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y/zW(l,l,l), u ( 2 n) 
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where 
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The average magnetization m for the nth generation 
is then given by 



straightforward calculations (see more detail [27j) show 
that one has 
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where 8 p — exp(/3J p ), #i = exp(/3Ji). 
We rewrite ([3]) in the reduced variables 
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Ml + 7i4 
M2 ~ U 5 

U± + Ui 



(3) 



where 7'"' = 2/1 + y 2 l> + 2/3" ; ■ It is quite obvious to 
note from (jH that the set D = {(a;, 0,0,0) : x E R + } is 
invariant with respect to that dynamical system. In this 
case the system is reduced to the following one: 



f(x) = 



1 
2B[ 



9 p + l)a; + 2) 



(6) 



Given the conditions P < 1 and 6± > 1 only one fixed 
point x* of f(x) can be found. 

The derived recursion relations ^ provide us (nu- 
merically) with the exact phase diagram in the 
(T/Ji, —Jp/Ji) space. Starting from random initial con- 
ditions (subject to the constraint j/1,2/2, 2/3 7^ 0); we mav 
observe the behavior of the recurrence relations ^ af- 
ter a large number of numerical iterations. The phases 
are characterized by the sequence of stable points of the 
recursion relations. Namely, in the simplest case a fixed 
point (x* , y{ , y 2 , y%) is reached. This point corresponds 
to a paramagnetic phase if y\ — y 2 = y 3 =0, or to a 
ferromagnetic phase if y\, y 2 ,y 3 =/= 0. The lower-order 
commensurate phases (short period) are described by 
a sequence of a few fixed points while the higher-order 
commensurate phases (large period) or incommensurate 
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FIG. 1: Phase Diagram of the model: F means ferromagnetic, 
P paramagnetic, M modulated, PM paramodulated, and < 
2 > the antiphase structure, i.e. the structure with period 4. 



phases (infinite period) are described by the quasicontin- 
uous or continuous attractors, respectively. The distinc- 
tion between a truly aperiodic case and one with a very 
long period is difficult to make numerically. 

The obtained phase diagram is presented in Fig. [1] 
At T — only two different ground-states are encoun- 
tered: the ferromagnetic states for —J p /Ji < 1/3, and 
states of period 4 (for example states with structure 
(1122), (1133), (2233)) for -J p /Ji > 1/3. Such states are 
called antiphase and denoted as < 2 >. The main fea- 
ture to be noted is the existence of a multiphase at finite 
temperature (T = 0, — J p /Ji = 1/3); this is where the 
paramagnetic (P), ferromagnetic (F), antiphase < 2 > 
(see below), modulated (M) (this phase contains both 
commensurate and incommensurate regions) and (PM) 
paramodulated phases, meet. We note the finding of a 
new phase wc refer to as paramodulated, which is char- 
acterized by 2-periodic points of the dynamical system 
© and lies inside the modulated phase. From |(SJ) one 
can see that in this region the average magnetization is 
the same as in the Paramagnetic one. It is important 
to note that such a phase is absent in the case of the 
Ising model, thus exhibiting a siginificant difference be- 
tween phase diagrams for Potts and Ising models. The 
rest of the diagram is quite similar to the one obtained 
by Vannimenus 1J] for the Ising model with similar in- 
teractions. 

Below we detail out the critical lines encountered in 
the phase diagram. 

The Paramagnetic Phase. The transition lines of the 
para-ferro and para-modulated are found to be continu- 
ous. Such lines are obtained by linearizing the system (j4|) 
around the fixed point (x*, 0, 0, 0). Note that the param- 
eters j/i, 2/2 and 2/3 vanish in the region R. The variable 



x* is unaffected in first order in 2/1, 2/2 and 2/3- 

The eigenvalue equation of the linearized system has 
the following form: 



A 3 + 
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The fixed point is linearly stable if the eigenvalues has 
moduli smaller than one. It is interesting to note that 
in the Paramagnetic case from ([5]), one may observe the 
average magnetization value as being equal to 2. 

To find out the transition lines it is necessary to exam- 
ine several cases with respect to whether the the eigen- 
values are real or complex. 

The Para-Ferro Transition. When the eigenvalues are 
real, the transition line will be characterized by the crite- 
rion that the largest (in absolute value) eigenvalue should 
be equal to unity. This determines the stability limit line 
we were looking for (since all fixed points are linearly 
stable if the eigenvalues have moduli smaller than one). 

The condition that the largest eigenvalue is equal to 1 
becomes 



(60 p + 2)x* 
~(49' 2 + 89 



' + (30| 
- 12)x* 



129 2 - 79 p + 28)x 
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- 46»p = 
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with x* positive. From this equation one can find x* and 
taking into account that is a fixed point of ([6]), one gets 
an equation of the para-ferro transition line which will 
be in the form 9\ = g(9 p ) for some function g. 

Note that the case 9 p = 1 (i.e. J p — 0) corresponds 
to the simple Potts model with nearest neighbor inter- 
actions and one recovers the well-known result for the 
critical temperature: exp(J/T c ) = 4 (jl. Hoi]. 

Observations of ^ show that at low temperatures 
{9 p < 1) from © one has x* ~ 9 p /3 and 9 X ~ (27/4)6»" 3 . 
In terms of T and p = —J p /J, the equation of the tran- 
sition line is given by 



27 

I -3p = Tlog ( — 



(T->0), 



which is in agreement with the slope obtained numeri- 
cally (see Fig. Ql. 

The Para- Modulated Transition. When the three 
eigenvalues are a pair of complex conjugates and one 
real, then the fixed point is approached in an oscilla- 
tory way and stability is achieved if the absolute values 
of all eigenvalues are less than 1. The critical transition 
line will then be characterized by the criterion that the 
modules of the complex eigenvalues are equal to unity. 
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Therefore, the instability occurs when eigenvalue A = i. 
Hence, the eigenvalue equation (JT]) is reduced to 



(56% + 29 p - 3)x* 3 + (56% + 6flg + 9 p )x 
+4(9 2 p + 39 p - l)x*+49 p = 0. 



(9) 



If 56*2 + 29 p - 3 > 0, i.e. 9 p > §, then equation (J2J 
does not allow for any positive solution. Therefore, the 

3 
5 : 



transition exists only if 9 P < | , that is -f- < log | , 



J P T 5 

> — log - 

J J 3 



and it corresponds to the asymptote of the transition line 
for large T in Fig. [T] 

Observations of @ show at low temperatures (9 p -C 
1), x* ~ d p and from © one finds 6 X ~ (l/2)6>- 3 . In 



terms of T and p 
line becomes 



-J p / J, the equation of the transition 



l-3p=--lo g 2 



(T-0) 



which is in agreement with the slope obtained numeri- 
cally. 

TTie Paramodulated Phase. 



between Potts and Ising models 14j, [l7[ • In Fig. [3] and 
Fig. [4] plots of the values of magnetization m are given at 
selected values —J p /J = 0.36 and — J p /J — 0.56, respec- 
tively. Note that this phase also issues from the multi- 
critical point (the Lifshitz point) T = 0, —J p /Ji = 1/3. 
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FIG. 3: The average magnetization versus T/J for — J p /J 
0.36. 




FIG. 2: Enlargement of Fig. [TJ around the PM phase 

As noted earlier that the set D = {(x, 0, 0, 0) : x G 
is invariant for the dynamical system one can con- 
sider its restriction to this set, which is given by ([5]). Nu- 
merical investigations of the dynamics of © show that 
in some values of the parameters T j J\ and — J p / J\ it has 
2-periodic points (note that no other periodic points are 
present). In the plane (T/Ji, J p /J\) such values form an 
island inside the modulated region (see Fig. [2]). In this 
island, the average magnetization m is equal to 2 as in the 
paramagnetic case. Therefore, we denote it as paramod- 
ulated (PM) phase. A similar phase is not present in 
the case of the Ising model, hence marking a difference 
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4: The average magnetization versus T/J for —J p / J 



Discussions. We have found an exact phase diagram of 
the Potts model with next nearest neighbor interactions 
on the Bethe lattice of order two. The diagram consists 
of five phases: ferromagnetic, paramagnetic, modulated, 
antiphase and paramodulated, all meeting at the Lifshift 
point i.e. p — 1/3. A distinctive feature of the diagram is 
seen in the existence of a new phase called the paramodu- 
lated phase found at low temperatures and characterized 
by 2-periodic points of the dynamical system ([6]). This 
phase lies inside the modulated phase and is inherent in 
the Potts model as no analogue can be found within the 
Ising setting. We believe the appearance of such a phase 
to essentially be another form of symmetry resulting in 
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the increasing of the number of spin from two to three. 
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